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HYBRID SIMULATION OF A CONTROL SYSTEM FOR A 
TUBULAR CHEMICAL REACTOR 



INTRODUCTION 



These Notes describe the application of hybrid techniques in the 
simulation of a distributed parameter system — specifically, the 
control system for a tubular chemical reactor—consisting of a 
number of partial differential equations. The simulation is per- 
formed on the EAI HYDAC® 2400 Hybrid Computing System which 
combines analog speed with digital memory and accuracy to pro- 
vide maximum problem solving capability of the available com- 
puting facilities. 

Figure 1 shows a block diagram of the physical system simulated. 
This system differs from other design, control or data analysis 
problems solved on analog computers alone by the chemical 
process industries (CPI) in that its mathematical description 
involves partial differential equations (PDE). Since each PDE 
must be transformed into a set of ordinary differential equations, 
each of which must be solved simultaneously, the implementation 
of a large number of almost identical circuits is required and this 
usually is sufficient to exhaust the analog equipment complement 
of even the best equipped laboratories. 

By means of a different mathematical formulation of the equations — 
a formulation which requires a function memory capability — this 
situation can be changed. Hybrid computers have satisfied this 
requirement for memory and thus enable system analysts to utilize 
the full power of available computing hardware. In addition, they 
have provided simple means for implementing logical operations, 
transport delays, and the transient simulation of stage-wise 
processes. 
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reactor, a separator, and a control loop. The 
kinetics of the system are: 



Figure 1: Block Diagram of Physical System 



OBJECTIVES 

The objectives of this simulation are typical of 
those encountered in many CPI simulations. These 
are: 



(1) To define the physical system mathe- 
matically, 



(2) To propose a control system to maintain 
a constant output product flow rate, and 



A + B C + D 

k„ 
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C + B ^ E + D 
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where component E, the desired product, is 
formed by feeding A, B and C to the reactor. 
Component D is a waste product which is dis- 
carded. 

Disturbances (step changes and noise such as 
valve chatter or low frequency sinusoids) enter 
the system through the make-up feed stream, 
A i? and cause variations in product flow, Ej. 
The purpose of the control loop is to operate on 
the flow error, e , using a control system to vary 
the feed rate of component C. Since B is present 
in excess (to maintain favorable reactor behavior) , 
variations in C, as in equation (1), have a direct 
effect on product formation. The control system 
selected is a conventional three-mode controller. 

It has been assumed that the separator is 100% 
efficient and that its dynamics can be represented 
by a first-order lag. Sensor and valve dynamics 
associated with the control loop are assumed 
negligible and the valve-position- versus-flow re- 
lationship is considered to be linear. 



MATHEMATICAL MODEL 



(3) To determine controller settings for the 
proposed control system. 



Kinetics: The rates of conversion for the com- 
ponents are given by: 



Problem definition, organization, programming, 
and hardware requirements to attain these objec- 
tives are presented in detail. 



r a = -k^ ab + k^ cd 



(3) 



r^ = -k.^ ab + k^ cd - k^ be + k^ de (4) 



PHYSICAL SYSTEM 



r^ = k x ab - k g cd - k 2 be + k 4 de (5) 



The complete system under study, shown sche- 
matically in Figure 1, consists of a tubular 



r d = k 1 ab-k 3 cd + k 2 bc-k 4de 



(6) 
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r = k rt be - de 
e 2 4 



(7) Control System: The classical transfer function 

for a three-mode controller is 



Where a...e are component concentrations, r a 
...r e are rates of formation, and k-^..^ are 
reaction rate constants. 



Reactor Equations: The concentration dependence 
of components A and E with respect to time and 
position is expressed by two partial differential 
equations: 



€ C l<xr D S + 1 




(16) 



where r D is the derivative time constant, r r is 
the reset time constant, is derivative gain, K r 
is reset gain, K c is controller gain, and A£is the 
change in percent of valve opening. 



dA T7 6a 
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de , y de 
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(9) 



where t is time, x is distance from reactor en- 
trance, and V is velocity of mixture through 
reactor. 



These equations can be rewritten in terms of flow 
rates to obtain 



Feed and Recycle Equation: Based on the linear valve 
assumption, the following reactor feed equations 
are obtained: 



A G= A 1+ A g+ f 
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The disturbance, f, can take two forms, viz: 



a step change, f =K 



(19) 



dA dA ^ + k — 

at ax IF 3 F 
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or a sinusoid, f = K sin cot 



(20) 
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B = B -A +A-E 
o o 



C=C + A - A - E 
o o 



D = B - B 

o 



(11) 

(12) 
(13) 
(14) 



Computational Model: Since the reactor- separator 
system is presumed to be operating at its design 
point when a disturbance is introduced into the 
system, good scaling dictates that equations be 
rewritten in perturbation form. The non-lineari- 
ties in the reactor equations make perturbation 
variables impractical. However, their form does 
lend itself to solution by the method of character- 
istics. 



Separator Equations: The separator equations, 
which are first-order lags, are, in transfer func- 
tion notation, 



Applying this technique yields 

dA 
"dy 



dA 

-- = kj_ AB - k CD 



(21) 
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E rS+1 
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(15) 



where S is the Laplace operator and r is the 
separator time constant. 



= k DE - k BC 
dy 4 2 



where 
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These reactor equations, solved in conjunction 
with the three algebraic equations, (12), (13), 
and (14), defining B, C and D in terms of A and 
E are the computational model of the reactor. 
The controller equation, (16), requires no modi- 
fication since it is already in perturbation form. 



A high-speed computational procedure is required 
to allow the reactor outlet product flows to be 
obtained in a semi -continuous manner. 

Please send for Application Study: 6.4.2h, Bul- 
letin No. ALHC 65076 for complete details, in- 
cluding results, of this simulation. 
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